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Abstract. Aims. In this paper we present the chemical abundance analysis from high resolution UVES spectra of seventeen 
bright giant stars of the Globular Cluster (GC) M 22. 

Results. We obtained an average iron abundance of [Fe/H] = -1.76 ± 0.02 (internal errors only) and an a enhancement of 
0.36 ± 0.04 (internal errors only). Na and O, and Al and O follow the well known anti-correlation found in many other GCs. 
We identified two groups of stars with significantly different abundances of the j-process elements Y, Zr and Ba. The relative 
numbers of the two group members are very similar to the ratio of the stars in the two SGBs of M 22 recently found by Piotto 
(2009). Y and Ba abundances do not correlate with Na, O and Al. The .s-element rich stars are also richer in iron and have 
higher Ca abundances. 

The results from high resolution spectra have been further confirmed by lower resolution GIRAFFE spectra of fourteen addi- 
tional M 22 stars. GIRAFFE spectra show also that the Eu - a pure r-process element - abundance is not related to the iron 
content. We discuss the chemical abundance pattern of M 22 stars in the context of the multiple stellar populations in GC 
scenario. 
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1. Introduction 

Globular Clusters (GC) are generally chemically homogeneous 
in their Fe-peak elements, while they show star-to-star abun- 
dance variations in light elements like C, N, O, Na, Mg, Al 
among others. In some cases, these chemical inhomogeneities 
result in well defined anticorrelations. For example, all GCs for 
which Na and O abundances have been measured, show a well 
defined NaO anti-correlation (see Carretta et al. 2006, 2008 for 
an update), often associated with an anti-correlation between 
Mg and Al contents. The origin of these variations has not yet 
been well understood, since both a primordial and an evolution- 
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ary explanation, or a combination of both have been proposed 
(see Gratton et al. 2004 for a recent review). 

Interestingly enough, abundance anomalies have been 
found also among stars in the lower part of the red giant 
branch (RGB) or even in the main sequence (MS). E.g., Cannon 
et al. (1998) found a CN bi-modality in MS and sub giant 
branch (SGB) stars in 47 Tuc, and the NaO anti-correlation 
was observed at the level of MS Turn Off (TO) and SGB in 
M 13 (Cohen & Melendez 2005), NGC 6397 and NGC 6752 
(Carretta et al. 2005; Gratton et al. 2001), NGC 6838 (Ramirez 
& Cohen 2002). By themselves, these results suggest the possi- 
bility of a primordial origin for the abundance inhomogeneities. 

More recent results, made possible by a significant im- 
provement of the photometric precision on HST images, in- 
dicating the presence of multimodal sequences in the color- 
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magnitude diagrams (CMD) of some GCs (Bedin et al. 2004, 
Piotto et al. 2007, Milone et al. 2008, and Piotto et al. in prepa- 
ration) further confirm that at least in some GCs there is more 
than one generation of stars formed from material chemically 
contaminated by previous generations. 

Among the clusters with multiple stellar populations, 
a) Centauri is the most complex and interesting case. This ob- 
ject is the only one for which variations in iron-peak elements 
have been certainly identified (Freeman & Rodgers 1975; 
Norris et al. 1996; Suntzeff & Kraft 1996). The Fe multimodal 
distribution is at least in part responsible of the multiple RGB 
(Lee et al. 1999, Pancino et al. 2000) of a> Cen. Also the MS 
of this GC splits into three sequences as shown by Bedin et al. 
(2004). Among the two principal MSs, the bluest one is more 
metal rich than the redder one (Piotto et al. 2005). So far, the 
only way we have to explain the photometric and spectroscopic 
properties of the MS of oj Cen is to assume that the bluest MS 
is also strongly He-enhanced. Recently, Villanova et al. (2007) 
showed that also the SGB splits in at least four branches with 
a large age difference (larger than 1 Gyr) among the different 
populations. 

NGC 2808 is the second (in time) cluster in which a MS 
split into three branches was found (Piotto et al. 2007). Also in 
this case, the MS multimodality was associated with different 
values of Helium, to the observed multimodal distribution of 
the stars along the NaO anti-correlation, where three groups of 
stars with different O (and Na) content were found by Carretta 
et al. (2006). 

In NGC 1851 Yong & Grundahl (2008) found abundance 
variations in various elements by studying a sample of 8 RGB 
stars. Sodium and oxygen follow the NaO anti-correlation. 
There is also some evidence for the presence of two groups of 
stars with different s-process element content, as well as of two 
groups of stars with different CN-strenght, that are possibly re- 
lated to the two sequences photometrically identified by Milone 
et al. (2008) along the SGB. The split SGB of NGC 1851 can 
be explained as due to the presence of two stellar populations. 
The stars in the fainter SGB could either be part of a popula- 
tion about 1 Gyr older than the bright SGB one, or could indeed 
be slightly younger than the bright SGB ones, but strongly en- 
riched in total C+N+O content (Cassisi et al. 2008). 

Another recent evidence for a primordial origin of abun- 
dance variations related to the presence of different populations 
of stars comes from Marino et al. (2008). By studying a large 
sample of RGB stars in the GC M 4, they found two distinct 
groups of stars with a different sodium content, which also dis- 
play a remarkable difference in the strenght of the CN-band. 
These two spectroscopic groups were found to populate two 
different regions along the RGB. They also noted that the RGB 
spread is present from the base of the RGB to the RGB-tip, 
suggesting that the spread must be related to the presence of 
two distinct stellar generations. 

At the basis of the present investigation, there is a very re- 
cent result by our group, who identified a bimodal distribution 
of the stars in the SGB of the GC NGC 6656 (Piotto 2009, see 
Fig.[l8]in the present paper). 

Located at a distance from the Sun of ~3.2 kpc (Harris 
1996), NGC 6656 (M 22) is a particularly interesting GC, be- 



cause a large number of photometric and spectroscopic studies 
suggested a complex metallicity spread, similar to, albeit sig- 
nificantly smaller than, that found in oj Centauri. In particular, 
it has been often suggested, though never convincingly con- 
firmed, that M 22 may have also a spread in the content of iron 
peak elements. 

The first evidence for a spread in metallicity comes from the 
significant spread along the RGB (Hesser et al. 1977; Peterson 
& Cudworth 1994) observed both in (B - V) and in Stromgren 
colors. However, it is still uncertain whether this spread can 
be attributed to a metallicity spread or to reddening variations. 
Due to its location on the sky, close to the Galactic plane and 
toward the Galactic Bulge [(b, €) (10°, 7°)], M 22 is affected 
by high and spatially varing interstellar absorption, with a red- 
dening in the interval 0.3< E(B - V) <0.5. This differential 
reddening creates a degeneracy in measuring metallicity when 
the atmospheric parameters of the stars are derived from their 
color. Spectroscopic studies are divided between those which 
conclude that no significant metallicity variations is present in 
M 22 (Cohen 1981, based on 3 stars; Gratton 1982, 4 stars) and 
studies claiming a spread in iron, with -1.4 < [Fe/H] < —1.9 
(Pilachowsky et al. 1984, 6 stars; Lehnert et al. 1991, 4 stars). 

Particularly interesting are the findings on CN-band 
strengths. Norris & Freeman (1983) showed that CN variations 
in M 22 were correlated with Ca, H and K line variations, sim- 
ilar to those in oj Cen. By studying a sample of 4 stars, Lehnert 
et al. (1991) found Ca and Fe variations that also correlated 
with variations in CH and CN band strength. However, Brown 
& Wallerstein (1992) found no Ca abundance differences be- 
tween CN-strong and CN-weak stars, though they observed 
differences in [Fe/H] correlating with the CN-strength. More 
recently, Kayser et al. (2008) found some indications of a CN- 
CH anti-correlation in SGB stars, maybe diluted by large un- 
certainties introduced by differential reddening. 

In the present study we analyze high resolution UVES 
spectra for a sample of seventeen RGB stars in M 22 in order 
to study the chemical abundances and possible relations with 
the recently found SGB split. In order to increase the statisti- 
cal significance and reinforce our findings, we also added the 
results from a sample of fourteen RGB stars from medium res- 
olution, high S/N GIRAFFE spectra. In Section [2] we provide 
an overview of the observations and of the data analysis, and in 
Section[3]we describe the procedure used to derive the chemical 
abundances. Our results on the chemical composition of M 22 
are presented in Section|U and a discussion on them is provided 
in Section[5] In Section|6]we look for possible connections be- 
tween our spectroscopic results and the two stellar populations 
photometrically observed by Piotto (2009). A comparison be- 
tween the results of this paper and those of Marino et al. (2008) 
on the GC M 4 is provided in Section|7] In Section[8]we present 
a brief discussion on the results obtained from GIRAFFE spec- 
tra. Section [9] summarizes the most relevant properties of the 
two stellar populations of M 22. 

2. Observations and membership analysis 

Our data set consists of spectra of seventeen RGB stars re- 
trieved from the ESO archive. The observations were ob- 
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ID RA [degree] DEC [degree] RV H [kms-'] B V I 



71 


279.032271 


-23.848980 


-136.74 


13.7458 


12.3171 


10.6334 


88 


279.150166 


-23.837640 


-151.21 


13.9042 


12.5003 


10.9469 


51 


279.107744 


-23.932690 


-166.65 


13.6157 


12.0413 


10.3165 


61 


279.087911 


-23.945560 


-145.81 


13.6269 


12.2122 


10.5598 


224 


279.056427 


-23.915140 


-147.87 


14.6701 


13.4704 


12.0023 


221 


279.137364 


-23.916380 


-156.45 


14.6746 


13.4606 


11.9562 


200043 


279.133958 


-23.934417 


-152.90 


13.4254 


11.9082 


10.1801 


200025 


279.042375 


-23.906056 


-153.22 


13.2488 


11.5217 


9.6761 


200068 


279.133875 


-23.858667 


-139.10 


13.7589 


12.2961 


10.6084 


200031 


279.113999 


-23.857222 


-142.08 


13.1504 


11.6354 


9.9282 


200076 


279.085417 


-23.940111 


-158.38 


13.7089 


12.3647 


10.7539 


200101 


279.159833 


-23.900361 


-131.00 


14.0444 


12.6726 


11.0746 


200080 


279.136750 


-23.852972 


-147.42 


13.7257 


12.4310 


10.8985 


200104 


279.083958 


-23.931611 


-133.25 


13.9028 


12.6644 


11.1595 


200083 


279.127417 


-23.880083 


-129.00 




12.4350 


10.8732 


200006 


279.072917 


-23.907249 


-148.16 


12.8797 


11.0304 


9.0975 


200005 


279.171042 


-23.971999 


-149.84 


12.7874 


10.9295 


8.9453 



tained using UVES (Dekker et al. 2000) and FLAMES @UVES 
(Pasquini et al. 2002). The spectra cover the wavelength range 
4800-6800 A, have a resolution R^45000, and have a typical 
S/N ~ 100 - 120. 

Data were reduced using UVES pipelines (Ballester 
et al. 2000), including bias subtraction, flat-field correc- 
tion, wavelength calibration, sky subtraction, and spectral 
rectification. The membership of the analyzed stars was 
established from the radial velocities obtained using the 
IRAF@FXCOR task, which cross-correlates the object spec- 
trum with a template. As template, we used a synthetic spec- 
trum obtained through the spectral synthesis code SPECTRUM 
(see http://www.phys.appstate.edu/spectrum/spectrum.html 
for more details), using a Kurucz model atmosphere with 
roughly the mean atmospheric parameters of our stars 
T eff = 4500 K, log(g) = 1.3, v t = 1.6 km/s, [Fe/H] = -1.70. At 
the end, each radial velocity was corrected to the heliocentric 
system. We obtained a mean radial velocity of -146 ± 2 km/s 
from all the selected spectra, which agrees well with the val- 
ues in the literature (Peterson & Cudworth 1994). Within 2cr, 
where cr is our measured velocity dispersion of 10 km/s, all our 
stars are members. The list of the analyzed stars, their coordi- 
nates, radial velocities, and magnitudes, are reported in Tabl^TJ 
Figure 1 shows the location of the target stars in the CMD of 
M22. 

We also analysed a sample of stars observed with 
GIRAFFE HR09, HR13, and HR15 set-ups at a resolution 
of R~20000-25000. These spectra were reduced by using the 
pipeline developed by Geneva observatory (Blecha et al. 2000). 
More datails on these data and their analysis are provided in 
Section[8] 

3. Abundance analysis 

Abundances for all elements, with the exception of oxygen, 
were measured from an equivalent width analysis by using the 
Local Thermodynamical Equilibrium (LTE) program MOOG 
(freely distributed by C. Sneden, University of Texas at Austin). 



The atmospheric parameters, i.e. temperature, gravity, and 
micro-turbolence, were determined from Fe lines by removing 
trends in the Excitation Potential (EP) and Equivalent Widths 
(EW) vs. abundance respectively, and satisfying the ionization 
equilibrium. 

At odd with other elements, we measured O content by 
comparing observed spectra with synthetic ones, because of the 
blending of the target O line at 6300 A with other spectral fea- 
tures. 

More details on the line-list, atmospheric parameters and 
abundance measurements can be found in Marino et al. (2008). 
In Table [2] we report the reference chemical abundances ob- 
tained for the Sun used in this paper. The obtained stellar pa- 
rameters for the M 22 analyzed stars are listed in Table[3] 



Table 2. Measured solar abundances 

(loge(X) = log(N x /N H ) + 12). 



Element 


abund. 


# lines 


Element 


abund. 


# lines 


O 


8.83 


1 


Tiiiii 


4.96 


12 


Na 


6.31 


4 


V 


3.89 


17 


Mg 


7.54 


3 


Fe FcI 


7.48 


145 


Al 


6.43 


2 


Fe FcI1 


7.51 


14 


Si 


7.61 


13 


Ni 


6.26 


47 


Ca 


6.39 


16 


Yyii 


2.24 


8 


SCsdl 


3.12 


12 


Zrzrii 


2.37 


1 


TlTil 


4.94 


33 


Basaii 


2.45 


2 



It is important to remark here that the method employed for 
the measurements of the atmospheric parameters is based on 
the spectra, and hence our temperatures are not color depen- 
dent. This is an important advantage in analyzing GC, as M 22, 
affected by high differential reddening. 

In Fig. [TJ we marked with red symbols the location of 
our seventeen UVES target stars on the B vs. (B - I) CMD. 
Photometry has been obtained with the Wide Field Imager 
(WFI) camera at the ESO/MPI 2.2 m telescope by Monaco 
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Table 3. Atmospheric parameters for the analyzed stars. 



ID 


Teff [K] 


log(g) 


v t [km s 1 ] 


[Fe/H] [dex] 


71 


4460 


1.15 


1.44 


-1.76 


88 


4460 


1.15 


1.55 


-1.70 


51 


4260 


0.90 


1.60 


-1.63 


61 


4430 


1.05 


1.70 


-1.78 


224 


4700 


1.70 


1.45 


-1.76 


221 


4750 


1.66 


1.20 


-1.75 


200043 


4400 


1.01 


1.70 


-1.77 


200025 


4100 


0.67 


1.80 


-1.62 


200068 


4500 


1.30 


1.52 


-1.84 


200031 


4300 


0.77 


1.55 


-1.85 


200076 


4500 


1.23 


1.35 


-1.83 


200101 


4500 


1.35 


1.55 


-1.74 


200080 


4600 


1.00 


1.45 


-1.81 


200104 


4700 


1.35 


1.75 


-1.92 


200083 


4490 


1.46 


1.66 


-1.63 


200006 


3990 


0.20 


2.08 


-1.66 


200005 


4000 


0.05 


2.02 


-1.94 



Table 4. The average abundance for M 22 stars. 



El. 


Abundance [dex] 


°" obs 


N star8 


[O/Fe] 


0.28±0.05 


0.20 


17 


[Na/Fe] 


0.24+0.08 


0.30 


17 


[Mg/Fe] 


0.39±0.03 


0.11 


17 


[Al/Fe] 


0.34+0.08 


0.31 


17 


[Si/Fe] 


0.43+0.01 


0.03 


17 


[Ca/Fe] 


0.31+0.02 


0.07 


17 


[Sc/Fe] ScII 


0.04+0.01 


0.04 


17 


[Ti/Feki 


0.24+0.01 


0.06 


17 


[Ti/Fe] T m 


0.34+0.01 


0.06 


17 


[V/Fe] 


-0.09+0.02 


0.10 


17 


[Cr/Fe] 


-0.13+0.02 


0.08 


17 


[Fe/H] 


-1.76+0.02 


0.10 


17 


[Ni/Fe] 


-0.07+0.01 


0.04 


17 


[Y/Fe] YII 


0.05+0.07 


0.27 


17 


[Zr/Fe]zrii 


0.36+0.06 


0.23 


17 


[Ba/Fe] BalI 


0.19+0.06 


0.23 


17 



et al. (2004) and stars with the highest photometric quality 
were carefully selected. Magnitudes have been corrected for 
sky concentration ( i.e. the increasing of the background level 
in frames near the center due to light reflected back from the de- 
tector to the optics.) by using the best available solution (Bellini 
et al. 2009). Since the photometry is affected by spatially vari- 
able interstellar reddening, we have corrected the CMD for 
this effect (by using the method described in Sarajedini et al. 
2007). To separate probable asymptotic giant branch (AGB) 
from RGB stars, we have drawn by hand the black dashed line 
of Fig. [TJ Stars that, on the basis of their position in the CMD, 
are possible AGB stars, will be marked from here on as trian- 
gles, while RGB stars with circles. 
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Fig. 1. Distribution of the UVES target stars (red symbols) on 
the B vs. (B - I) CMD corrected for differential reddening. 
The star #200083 is not plotted because the B magnitude is 
not available for this star. The stars observed with GIRAFFE 
are shown by blue crosses. 

3.1. Internal errors associated to chemical 
abundances 

The main goal of this paper is to study the intrinsic variation of 
chemical abundances affecting M 22 stars. 

The differences in the measured chemical abundances from 
star-to-star are a consequence of both measurements errors and 
intrinsic variations in their chemical composition. In this sec- 
tion, our final goal is to disentangle between internal errors, 
due to measurements uncertainties, and real variations in the 
chemical content of the stars. To this aim, we will compare the 
observed dispersion in the chemical abundances cr b s , listed in 
Table @] with that produced by internal errors cr tot . Since we 
are interested in the study of star-to-star intrinsic abundance 
variations, we do not treat possible external sources of errors 
that do not affect relative abundances. 

Two sources of errors mainly contribute to cr tot : the uncer- 
tainties in the EWs, and the uncertainties in the atmospheric 
parameters. 

In order to derive the typical error in the EWs, we consider 
two stars (#200101 and #71) with similar atmospheric parame- 
ters and roughly the same iron abundance. In this way, any dif- 
ference in the EWs in the iron spectral lines, can be attributed 
to measurement errors. The dispersion of the distribution of the 
differences between the EWs of the iron lines of the two se- 
lected stars, that is 2.3 mA, has been taken as our estimate of 
the typical error in the EWs measurement. The corresponding 
error in the chemical abundances has been calculated by vary- 
ing the EWs of a star (#200101) at intermediate temperature, 
representative of our sample, by 2.3 mA. The variations in the 
obtained abundances for each chemical specie, listed in Table[5] 
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(column 7), have been taken as our best estimate of the internal 
errors introduced by uncertainties in the EWs. 

We used the same procedure described in Marino et al. 
(2008) in order to estimate the uncertainties associated to the 
atmospheric parameters, and the corresponding errors related 
to chemical abundances. From our analysis, the obtained un- 
certainties related to atmospheric parameters are: AT e ff=+50 
K, Alog(g)=+0.14, and Av t =+0.13 km/s. These internal errors 
in the atmospheric parameters translate into the errors in chem- 
ical abundances listed in Table [5] (columns 2, 3 and 4). 

We investigated also the influence of a variation in the to- 
tal metallicity ([A/H]) of the model atmosphere on the derived 
abundances. By varying the metallicity of the model by 0.10 
dex, that is the iron observed dispersion, the element abun- 
dances change by the amount listed in Col. 5 of Table 5. A 
variation in the metallicity of the model atmosphere mainly 
changes the ionization equilibrium, and hence the values of 
log(g), since we used the ionization equilibrium between Fel 
and Fell to derive gravities. We calculated that by increasing 
[A/H] by 0.10 dex, log(g) decreases by ~0.06, while temper- 
ature and microturbolence do not change significantly. Since 
we are interested in the search for possible small star-to-star 
variations of iron abundances, we measured the effect in the 
[Fe/H] abundances due to this change in gravity in a model at- 
mosphere with increased metallicity, and verified that it does 
not affect the derived iron abundances by more than 0.01 dex. 
By increasing the total metallicity by 0.2 dex, the Fell abun- 
dances change by ~0.06 dex, and we have to decrease log(g) 
by ~0. 12 to re-establish the ionization equilibrium between Fel 
and Fell. Also in this case the effect on the iron abundances is 
smaller than 0.01 dex. 

Column 9 of Table 5 reports the quadratic sum (cr tot ) of the 
errors coming from the EWs (<x EW ) and from the atmospheric 
parameters (cr atm ) uncertainties. Column 8 gives the observed 
dispersion (cr obs ). 

Since the oxygen abundance was calculated from the same 
spectral line and the spectra have similar S/N, we assume as 
an estimate of the error related to [O/Fe], the <x tot calculated 
in Marino et al. (2008) for M 4 red giants. This error, for M 4 
stars, was calculated as the dispersion in O of the O-rich stars, 
i.e. the Na-poor group, assumed to be homogeneous in oxygen 
content. 

4. The chemical composition of M 22 

4.1. Iron-peak and a elements 

The wide spectral range of UVES data allows us to obtain 
chemical abundances for fifteen chemical species. Table|4]gives 
the mean abundance for each element (Col. 2), the rms of the 
mean of the abundances cr Q b s (Col. 3), and the number of stars 
(N s tars) used to calculate the mean (Col. 4). In Table |4] to each 
average abundance we associated an uncertainty which is the 
rms scatter (cr obs ) divided by VN s tars ~~ 1> although some of the 
distributions are clearly not Gaussians. A plot of our measured 
abundances is shown in Fig. |2l where, for each box, the cen- 
tral horizontal line is the mean value for each element, and 
the upper and lower lines contain the 68.27% of the distribu- 



tion around the mean. The points represent individual measure- 
ments. 

Our results confirm M 22 to be a metal poor GC with a 
mean metallicity: 

[Fe/H] =-1.76 + 0.02 dex, cr obs = 0.10 

where <x obs is the rms of the 17 measurements. The iron-peak 
elements Ni, Cr, and V have abundances of [Ni/Fe]=-0.07 + 
0.01, [Cr/Fe]=-0.13 ± 0.02, and [V/Fe]=-0.09 ± 0.02, respec- 
tively. 

We measured the chemical abundances for five a elements: 
O, Mg, Si, Ca and Ti. The corresponding abundances are listed 
in Table H] Here we consider only Mg, Si, Ca, and Ti. The re- 
sults for oxygen will be discussed in Section l4~2l These four a 
elements are all overabundant with respect to solar values, with 
an average of: 

[a/Fe]=+0.36 + 0.04 

For calcium we obtained a mean value of [Ca/Fe]=+0.31 ± 
0.02, similar to that found in other GCs, Interestingly enough, 
our stars show quite a large dispersion (cr obs = 0.07, see 
Table |4]i in the Ca abundance. This spread will be discussed 
more in detail in Section[5] 
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Fig. 2. Box plot for M 22 star elements abundances. The points 
are the individual measurements. The central horizontal line of 
each box is the mean of the data; the upper and lower lines 
contain the 68.27% of the distribution around the mean value. 



4.2. NaO anti-correlation 

As displayed in Fig. [3] sodium and oxygen show the typical 
NaO anti-correlation found in RGB stars in all GCs where Na 
and O have been measured so far (see Carretta et al. 2006). 
The [Na/Fe] values range from ~ -0.25 to ~0.7 dex, with a 
dispersion <x obs =0.30, and [O/Fe] abundances cover the interval 

from 0.10 to — 1-0.5 dex, with a dispersion <x obs =0.20. 

As it will be discussed in Section [6] Piotto (2009) have 
shown that the SGB of M 22 is split into two separate branches, 
indicating the presence of two stellar populations. In Fig. [4] we 
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Table 5. Sensitivity of derived UVES abundances to the atmospheric parameters and EWs. We reported the error cr atm due to the 
uncertainties in the atmospheric parameters (AT e( f, Alog(g), Av t , and A([A/H])) due to the error in EW measurements (<x EW ), the 
squared sum of the two (<x tot ), and the observed dispersion (cr obs ) for each element. 





AT ef) [K] 


Alog(g) 


Av, [km s '] 


A([A/H]) [dex] 


& atm 


C EW 


0"tot 


0"obs 




+50 


+0.14 


+0.13 


+0.10 










[O/Fe] 


- 


- 


- 


- 


- 


- 


0.04 


0.20 


[Na/Fe] 


-0.02 


-0.01 


+0.01 


-0.01 


0.03 


0.01 


0.03 


0.30 


[Mg/Fe] 


-0.03 


+0.00 


+0.01 


-0.01 


0.03 


0.02 


0.04 


0.11 


[Al/Fe] 


-0.03 


-0.01 


+0.02 


-0.01 


0.04 


0.06 


0.07 


0.31 


[Si/Fe] 


-0.05 


+0.01 


+0.01 


+0.01 


0.05 


0.01 


0.05 


0.03 


[Ca/Fe] 


-0.01 


-0.01 


-0.01 


-0.01 


0.02 


0.01 


0.02 


0.07 


[Sc/Fe] 


+0.03 


+0.00 


+0.00 


+0.00 


0.03 


0.03 


0.04 


0.04 


[Ti/Fe] ra 


+0.04 


-0.01 


+0.00 


-0.02 


0.05 


0.01 


0.05 


0.06 


[Ti/Fe] Ti ii 


+0.02 


-0.01 


-0.02 


-0.01 


0.03 


0.02 


0.04 


0.06 


[V/Fe] 


+0.03 


+0.00 


+0.03 


-0.01 


0.04 


0.01 


0.04 


0.10 


[Cr/Fe] 


+0.02 


-0.03 


-0.01 


-0.02 


0.04 


0.01 


0.04 


0.08 


[Fe/H] Fd 


+0.07 


+0.00 


-0.02 


-0.01 


0.07 


0.06 


0.09 


0.10 


[Fe/H] Fen 


-0.02 
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0.06 


0.06 


0.08 


0.10 


[Ni/Fe] 


-0.01 


+0.00 


+0.01 


+0.00 


0.01 


0.01 


0.01 


0.04 


[Y/Fe] 


+0.03 


-0.01 


-0.02 


+0.00 


0.04 


0.01 


0.04 


0.27 


[Zr/Fe] 


+0.03 


+0.00 


+0.02 


+0.00 


0.04 


0.01 


0.04 


0.23 


[Ba/Fe] 


+0.05 


-0.01 


-0.09 
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0.03 
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Fig. 3. The anti-correlation between [Na/Fe] and [O/Fe] abun- 
dance ratios. Triangles indicate stars that, on the basis of their 
position on the CMD, should be AGB. 

compare the Na versus O trend for M 22 with those observed 
in other four GCs where multiple stellar populations have 
been identified: NGC 2808 (data from Carretta et al. 2006), 
NGC 1851 (data from Yong & Grundahl 2008), NGC 6388 
(data from Carretta et al. 2007) and M 4 (data from Marino et 
al. 2008). 

In NGC 2808 the presence of multiple stellar populations 
is inferred by the triple MS (Piotto et al. 2007). The three MS 
branches may be associated with the complexities of the cluster 
horizontal branch (D' Antona & Caloi 2004) and of its oxygen 



abundance distribution (Carretta et al. 2006). On the contrary, 
the MS of M 22 is narrow and a spread or split, if any, must 
be smaller than 0.02 magnitudes in ;«F435w - WF814W color (see 
Piotto 2009). Figure|]shows that NGC 2808 stars cover almost 
the same range of Na abundances as M 22, but span a range of 
O abundance at least two times larger. 

In M 4, the presence of multiple populations is inferred 
by the bimodal distribution of the Na abundance (Marino et 
al. 2008). In addition, stars from the two groups with different 
Na content populate two distinct RGB sequences in the U vs. 
U-B CMD. Moreover, no MS or SGB split has been identi- 
fied. Interestingly, in M 4, the maximum variations of [Na/Fe] 
and [O/Fe] are respectively 0.39 and 0.31 dex smaller than 
those observed in M 22. In particular, a few O-poor stars with 
[O/Fe]<~0.00 and [Na/Fe] >~0.50 are present among M 22 
stars and lack in M 4. 

As in M 22, in NGC 1851 (Milone et al. 2008) and 
NGC 6388 (Piotto 2008; Moretti et al. 2008; Piotto 2009) 
the presence of multiple populations is inferred by a split of 
the SGB. Unfortunately, for both these clusters, the available 
chemical measurements from high resolution spectroscopy are 
limited to seven RGB stars for NGC 6388 (Carretta et al. 2007) 
and eight RGB stars for NGC 1851 (Yong & Grundahl 2008). 
In the case of NGC 6388, the stars are located in a portion of the 
NaO plane that is not populated by any M 22 star of our sample, 
being NGC 6388 stars systematically O-poorer. On the con- 
trary, the range of NaO anti-correlation in NGC 1851 matches 
quite well that of M 22. 

At variance with the case of M 4 and NGC 2808, both 
NGC 1851 and M 22 seem to show a continuum distribution 
in [Na/Fe], without hints of multi-modalities, also if the small 
number of stars studied in this paper and in Yong & Grundahl 
(2008) for NGC 1851, prevents us from definitively excluding 
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the presence of discontinuities in the [Na/Fe] or [O/Fe] distri- 
bution. 
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Fig. 4. NaO anti-correlation for M 22 stars superimposed to a 
collection of stars of four GCs that host multiple stellar popu- 
lations. 



4.3. Aluminum and Magnesium 

The abundances of aluminum and magnesium have been de- 
termined from Al lines at 6696 A, and 6698 A, from the Mg 
doublet at 63 1 8-63 1 9 A and the Mg line at 57 1 1 A. 

Figure [5] shows the [Mg/Fe] ratio as a function of [Al/Fe]. 
There is no clear MgAl correlation, despite of the presence of 
a well defined NaO anti-correlation, and of a clear AINa cor- 
relation, as shown in Fig. [6] Assuming that Na enhancement 
comes from proton capture process at the expenses of Ne, we 
expect to observe also a MgAl anti-correlation, since Al forms 
at the expenses of Mg. This means we expect a decrease of Mg 
abundance with the increasing of Al content. We do not ob- 
serve such an effect but, given our uncertainties, it could be too 
small to be detected. The lack of such a clear correlation was 
observed also in M 4 by Marino et al. (2008). However, they 
found a small difference in Mg content among stars character- 
ized by large Na content differences, according to the scenario 
proposed by Ivans et al. (1999), who predicted that a drop of 
only 0.05 dex in Mg is needed to account for the increase in 
abundance of Al (see their Section 4.2.2). 

4.4. s-process elements 

We have measured abundance for three s-process elements: yt- 
trium, zirconium, and barium. All of them span a wide range of 
abundance values. The maximum variations of [Y/Fe], [Zr/Fe] 
and [Ba/Fe] have amplitudes of 0.87, 0.73 and 0.65 dex respec- 
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Fig. 5. [Mg/Fe] vs. [Al/Fe] abundance ratios. 
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tively, despite of the small estimated internal error (cr tot < 0.1, 
see Table |3)- 

In Fig. [7] we show [Zr/Fe] and [Ba/Fe] as a function of 
[Y/Fe] (left and central panels), and [Ba/Fe] as a function of 
[Zr/Fe] (right panel). A clear correlation between each pair of 
i-process elements is evident, and also the slopes of the best 
fitting lines are similar. 

Most importantly, the .s-process elements clearly show 
a bimodal distribution: one group is overabundant in 
i-elements with average values of [Y/Fe]=+0.34+0.05, 
[Zr/Fe] =+0.60±0.04 and [Ba/Fe] =+0.43+0.04. This group 
contains seven out of seventeen stars (i.e. ~40% of the 
whole sample studied in this paper). The remaining ten 
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stars have [Y/Fe]=-0.16±0.03, [Zr/Fe] =+0.20±0.03, and 
[Ba/Fe] =+0.02+0.04. 

Figure [7] suggests that we can isolate a .s-process elements 
rich group, and a .s-process elements poor one, by selecting 
stars with [Y/Fe] greater and smaller than 0, respectively. Stars 
rich in .s-process elements are represented by red filled sym- 
bols, while s poor stars by the black empty ones. Red and 
black crosses with error bars indicate the average abundances 
for stars of each group. In Table [6] we have listed the mean 
abundances for all the elements studied in this paper, calcu- 
lated separately for each of these two groups. There are few 
remarkable differences in the average abundances for the two 
groups. 



Table 6. Average abundances for s-element process poor and 
rich stars. 



Element 


■s-process elements poor 


s-process elements rich 


[O/Fe] 


0.30±0.06 


0.26+0.10 


[Na/Fe] 


0.13+0.11 


0.40+0.08 


[Mg/Fe] 


0.36+0.04 


0.45+0.03 


[Al/Fe] 


0.26+0.11 


0.46+0.10 


[Si/Fe] 


0.42+0.01 


0.44+0.01 


[Ca/Fe] 


0.27+0.01 


0.38+0.02 


[Sc/Fe] 


0.04+0.02 


0.04+0.02 


[Ti/Fe] TU 


0.22+0.01 


0.27+0.03 


[Ti/FeJra, 


0.33+0.02 


0.36+0.03 


[V/Fe] 


-0.10+0.04 


-0.09+0.03 


[Cr/Fe] 


-0.16+0.02 


-0.09+0.04 


[Fe/H] 


-1.82+0.02 


-1.68+0.02 


[Ni/Fe] 


-0.09+0.01 


-0.06+0.01 


[Y/Fe] 


-0.16+0.03 


0.34+0.05 


[Zr/Fe] 


0.20+0.03 


0.60+0.04 


[Ba/Fe] 


0.02+0.04 


0.43+0.04 



As an example, in Fig. [8] we show the [Y/Fe] ratio as a 
function of the [Na/Fe], [Al/Fe] and [O/Fe] ratios. There is no 
clear correlation between Y, as well as Ba, and Zr abundances 
and Na, O or Al. However, stars enriched in .s-process elements 
are all Na-rich and Al-rich, while the group of stars with lower 



.s-process elements abundance spans over almost all the values 
of Na and Al abundance. 

The bi-modality in .s-process elements in M 22 resembles 
the case of NGC 1 85 1 . In NGC 1851, Yong & Grundahl (2008) 
noted that the abundances of the .s-process elements Zr and La 
appear to cluster around two distinct values. They suggested 
that the two corresponding groups of stars should be related to 
the two stellar populations photometrically observed by Milone 
et al. (2008) along the SGB. In NGC 1851 the s-element abun- 
dance appears also to correlate with Na, Al and O abundance. 




5860 



A[A] 

Fig. 9. The upper panel shows the spectra of the metal poor 
stars #200068 ([Fe/H] =-1.84) and of the metal rich star 
#200083 ([Fe/H] =-1.63) centered on the lines of Fel. The stel- 
lar parameters are very similar, however, the Fe lines differ 
significantly. The bottom panel shows the same lines for stars 
#200101 and #71 with similar metallicity and stellar parame- 
ters. 

By the similarity between our results and the ones by Yong 
& Grundahl (2008) on NGC 1851, it is tempting to associate 
the presence of the two groups of stars with different .s-process 
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Fig. 11. [Y/Fe] (left), [Zr/Fe] (center) and [Ba/Fe] (right) as a function of [Fe/H]. Symbols are as in Fig. [7] dashed lines are the 
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element content, with the two populations of stars isolated 
along the SGB of M 22 by Piotto (2009, see also FigfT8ll. 
Note that barium, yttrium, and zirconium can be considered 
as the signature of the s-process that occurs in intermediate 
mass AGB stars (Busso et al. 2001), whose wind could have 
polluted the primordial material from which the second gener- 
ation of stars in M 22 and NGC 1851 formed. Also the fact that 
i-element rich stars seem to be also rich in Na and Al could 



be due to the pollution from intermediate mass AGB stars of 
the material from which a second generation of stars formed, 
though the lack of a clear correlation is more difficult to inter- 
pret. 

We note that also the results by Piotto et al. (2005) on the 
metal content of the stars in the two main MSs of a> Centauri 
can be interpreted within this scenario. Piotto et al. (2005) 
found [Ba/Fe] -0.5 for the metal-poorer ([Fe/H] =-1.68) red 
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MS, and [Ba/Fe]~0.8 for the metal-richer ([Fe/H]=-1.37) 
bluer MS. This result was further confirmed by Villanova et 
al. (2007), who observed that SGB metal-poor stars have a Ba 
content lower than intermediate-metallicity ones by about 0.2 
dex. 

The relations between the ^-process element abundance and 
the iron content in M 22 is the argument of the following sec- 
tion. 

5. The spread in Fe of M 22 

As discussed in SectionQ] for long time, the existence of an in- 
trinsic Fe spread in M 22 has been debated in the literature (see 
Ivans et al. 2004 for a review), since photometric and spectro- 
scopic studies have yielded conflicting results. Some spectro- 
scopic studies found no significant variations (Gratton 1982, 
Ivans et al. 2004), whereas others seem to find a variation in 
[Fe/H] up to ~0.5 dex (Pilachowsky et al. 1984). Photometric 
studies gave similar controversial results. Undoubtedly, the 
RGB of M 22 has a large spread in color. This spread is ob- 
served both in BVI and Stromgren photometry, but interpreted 
either in terms of metallicity variation or differential reddening, 
or a combination of the two. 

We want to emphasize again that the [Fe/H] measurements 
presented in this paper are not based on photometric data, and 
therefore do not suffer the effects of differential reddening. For 
this reason, they constitute an appropriate tool to make the issue 
of [Fe/H] variations in M 22 clearer. 

From Table [5] we see that the observed dispersion <x obs in 
iron is comparable with the estimated internal error <x tot , i.e. the 
observed star-to-star metallicity scatter could be interpreted as 
due to measurement errors only. This fact demonstrates how 
difficult it is to establish the statistical significance of any in- 
trinsic spread in [Fe/H]. In this paper, we can tackle the prob- 
lem of the iron dispersion in M 22 in a different way. The ob- 
served dispersion in iron content alone could be a poor indica- 
tor of any intrinsic metallicity dispersion. First of all, because 
of the abundance measurements errors. In addition, we note 
that if anomalies in [Fe/H] affect only a small fraction of M 22 
stars, their effect on the iron dispersion of the whole sample of 
stars could result to be negligible. 

A visual inspection of some spectra reinforces the sugges- 
tion that there may be star-to-star iron variations. As an exam- 
ple, in Fig. [9] we show spectra of two pairs of stars with very 
similar stellar atmospheric parameters. The two spectral lines 
in Fig. [9] are iron lines. The upper panel shows the spectra of 
the stars #200068 and #200083 for which we have measured an 
iron abundance [Fe/H] =-1.84 and [Fe/H] =-1.63, respectively 
(see Table [3). The line depths differ significantly and, because 
of the similarity in the atmospheric parameter, must indicate an 
intrinsic iron difference. As a comparison, in the bottom panel 
we plot the same spectral region for two stars (#200101 and 
#71) with almost the same iron content. 

Our results on the iron dispersion can not be conclusive by 
themselves. However, it is rather instructive to look for possi- 
ble correlations between [Fe/H] and the other chemical abun- 
dances (mainly the ones showing a spread or bimodal distribu- 
tion, as the ^-process elements). 



Figure[lO]shows that Fe abundance is not significantly cor- 
related with the Al, Na, and O abundances, elements which, on 
the other hand, are involved in a well defined anti-correlation, 
as discussed in Section l4~2l 

However, when we compare the iron abundances with 
the ^-process element abundances, we find a strong corre- 
lation, with i-process element rich stars having systemati- 
cally higher [Fe/H], as shown in the three panels of Fig. QT| 
The two 5-process element rich and poor groups have av- 
erage [Fe/H\s-rich=- 1.68+0.02 and [Fe/H],_p OOr =-1.82±0.02, 
respectively (6 [Fe/H] =0.14+0.03). 

We want to emphasize that the significance in the Fe vari- 
ation can be appreciated only when we consider the average 
iron content of the two groups characterized by a different 
s-elements content. This is further demonstrated by a simple 
test we run, and whose results are summarized in Fig. [12] We 
simulated 100,000 stars with iron content, and iron abundance 
dispersion as the observed stars, i. e. 41% of the sample are 
i-process element rich stars (red Gaussian) and a 59% are s- 
process element poor ones (blue Gaussian) with a mean [Fe/H] 
of -1.82 and -1.68 dex, respectively. The dispersion in iron 
content of each group was taken equal to 0.09 dex, that is the 
error we estimated for the iron abundance (see Table 5). The 
resulting metallicity abundance dispersion for the total sample 
is of 0. 1 1 dex, which is close to, and only marginally greater, 
than the dispersion expected from the measurement errors. 




-2.2 -2 -1.8 -1.6 -1.4 

[Fe/H] 



Fig. 12. Simulation of 100,000 stars with the properties of the 
two observed groups of stars. The red and the blue histograms 
represent the s-rich and the s-poor groups of stars respectively. 

The presence of two groups of stars, one of which enriched 
both in i-process elements and iron, is analogous to the case of 
u> Centauri discussed in the previous Section |4~4l where Piotto 
et al. (2005) and Villanova et al. (2007) found that the more 
metal rich, He-enriched stars have also higher Ba content. 
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Fig. 14. [Y/Fe] {left), [Zr/Fe] {center) and [Ba/Fe] {right) as a function of [Ca/Fe]. Symbols are as in Fig. [7] dashed lines are 
the best least square fitting straight lines. 



0.5 



0.4 



0) 



ca 0.3 
o 



0.2 



0.1 



i i | i i i | i i i |- 

slope = 0.48 ±0.14 



o 



i -L i 
HjH, 



O 
© , 



A 



AS 



9- 
..T.cr 

o 

A 



1.8 



1.6 



[Fe/H] 
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As shown in Fig. Qj] the iron abundances are well corre- 
lated also with calcium abundances, and the Ca content cor- 
relates with the ^-process element abundance. Again, a sim- 
ilar behavior is present in a> Centauri, where a spread in Ca 
was known since Freeman & Rodgers (1975). More recently, 
Villanova et al. (2007) showed that the SGB population with 

[Fe/H] 1.2 has a mean [Ca/Fe] larger by ~ 0.1 dex, than that 

found in the more metal-poor population. 

In M 22, evidence for a calcium spread was already noted 
and Norris & Freeman (1983) showed a correlation between 
CN variations and Ca, similar to those in u> Centauri. Lehnert 
et al. (1991), by studying a sample of 4 stars, found both Ca 
and Fe variations correlated with the CN-band strengths. 

Our abundance measurements show that both calcium and 
iron correlate with the s-process elements (Fig.[14]and Fig.fTTTi. 
Moreover, as shown in Fig. Q3] we found that calcium, like 



iron, is not clearly correlated with Na, although the stars with 
an higher Ca content seem to be slightly Na rich. 

In Fig. 7, Fig. 11, Fig. 13, and Fig 14 six out of seven 
probable AGB stars, represented by triangles, belong to the s- 
poor group. We want to emphasize here that our selection of 
the probable AGB stars is based only on a visual inspection 
of the stars on the CMD, without considering photometric er- 
rors. In any case, assuming that all these stars are indeed real 
AGB members, our s-poor sample would include both RGB 
and AGB stars. If in our s-poor sample there is really a group 
of AGB stars, they could be the evolution of low mass stars of 
the primordial population, not able to activate the third dredge- 
up and enrich their surface of s-process elements. Hence they 
should trace the primordial composition of the cluster. Possibly, 
due to the larger iron content, stars with higher s element con- 
tent, have systematically redder colors with respect s-process 
element poor stars, and have low probability to be pushed by 
photometric errors in to the AGB region. 

Figure [16] shows [Mg/Fe] and [Si/Fe] as a function of 
[Fe/H]. Both Mg and Si are slightly overabundant in s-process 
element rich stars with respect to s-process element poor ones 
(see also Table [6]). This seems to suggest that core collapse 
SNe (CCSNe) are the best candidates to produce the iron ex- 
cess of the second generation of stars. Indeed, would iron be 
produced by Type la SNe, we would expect a lower [Mg/Fe] 
and [Si/Fe] ratio for the group of stars enriched in s-process 
elements with respect to the first generation of (s-process ele- 
ment poor) stars. In fact, SNela events are selectively enriched 
in iron (although modest quantity of Si are also produced). On 
the contrary, CCSNe along with iron, produce also Mg and Si 
in higher quantity with respect to SNela. 

With a present day mass of ~ 5 x 10 M (Pryor & Meylan 
1993) NGC 6656 is one of the most massive GCs of the Milky 
Way. The s-elements rich population, with a mass of ~ 2 x 10 5 
M and an iron abundance [Fe/H] =-1.68 ± 0.02 dex, includes 
~1.5 M Q of fresh iron if we assume Zg e = 0.0013. On average 
each CCSN produces -0.07 M Q of iron (Hamuy 2003), there- 
fore about twenty SNe are needed to produce the fresh iron of 
the second stellar population. In this scenario, the fainter SGB 
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Fig. 16. [Mg/Fe] (left), and [Si/Fe] (right) as a function of [Fe/H]. Symbols are as in Fig. [7] 
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Fig. 17. [Ba/Fe] (Ze/r), and [Ca/Fe] as a function [Fe/H] (rig/zf) for 104 stars in the GC M 4 (from Marino et al. 2008). 

6. Can the iron spread account for the SGB split? 
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Fig. 15. [Na/Fe] as a function of [Ca/Fe]. Symbols are as in 
Fig. E] 



(and TO) of this second generation of stars is attributed to their 
different chemical mixture, rather than to an age difference. 



Figure [18] from Piotto (2009) shows that the SGB of M 22 is 
splitted into two branches. In the previous section, we have 
shown that in M 22 there are two groups of stars, with differ- 
ent ^-process element contents, and with two different average 
iron contents. In this section,we want to investigate whether the 
different iron content can explain the split of the SGB. 

In Fig.|T9j we compare two isochrones from Pietrinferni et 
al. (2004) in the ACS/WFC plane w F 606w vs. m F 606w - ?«F8i4W 
of Fig. Q~8] Both of them have an age of 14 Gyr, but different 
metallicities. The black line corresponds to the mean metallic - 
ity of the group of ^-process element poor stars ([Fe/H]=- 1 .82) 
and the dashed red line is an isochrone with the average 
[Fe/H]=-1.68 of the s-process element rich stars. The differ- 
ent metallicity mainly reflects in a split of the RGB and of the 
SGB. In the inset, we show a zoom of the SGB region. At 
'«F606W - »JF8i4w = 0.85 the difference in magnitude ;«F606W 
between the two isochrones is £»jf606w = 0.10, about 0.07 
magnitudes smaller than that observed by Piotto (2009). We 
conclude that the observed difference in [Fe/H] can contribute 
to produce the splitting of the SGB observed by Piotto (2009), 
but it is not sufficient. On the other hand, the entire shape of the 
turn off-SGB-RGB region is difficult to reproduce with stan- 
dard, alpha-enhanced isochrones. Likely, this is due to the fact 
that the origin of the split may be much more complicated and 
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involves also the NaCNO abundances, as suggested by Cassisi 
et al. (2008) for the case of NGC 1851. 

7. Comparison with M 4 

In this section, we present a comparative analysis between the 
chemical abundances obtained in this paper for M 22 and the 
abundance measurements on M 4 RGB stars by Marino et al. 
(2008) to better outline the complexity of the multiple popula- 
tion appearence in different clusters. 

Such a comparison is rather instructive. First of all, because 
M 4, similarly to M 22, is affected by high differential redden- 
ing (Lyons et al. 1995, Ivans et al. 1999), and moreover, their 
spectra were analysed employing the same procedure. We note 
also that the UVES spectra of M 4 were collected with the same 
set-up and have almost the same S/N ratio as the M 22 spectra 
analysed in this work. 

Marino et al. (2008) have shown that M 4 hosts two distinct 
stellar populations, characterized by different Na content, and 
different CN-band strenght. These two groups of stars also de- 
fine two sequences along the RGB, but there is no SGB split. 
At variance with M 22, M 4 does not show any evidence of 
intrinsic Fe spread. Marino et al. (2008) set an upper limit for 
the [Fe/H] spread of 0.05 dex (lcr) in M 4. In the case of M 22, 
as discussed in Section 14.21 we identified a well defined NaO 
anticorrelation, but we have no evidence of a dichotomy in Na 
distribution, as in M 4. Instead, we found a dichotomy in the 
i-process elements, and the SGB is splitted into two branches. 

As a comparison, we show in Fig[l7]the Fe abundances as 
a function of Ba and Ca in M 4, using the same scale used for 
the same plot for our M 22 targets (Fig. QT| right panel, and 
Fig. [T3l respectively). In the case of M 4 there is no evidence 
of a chemical spread in [Ba/Fe] vs. [Fe/H] nor in [Ca/Fe] vs. 
[Fe/H] . Note that the rms of the calcium abundance in M 4 is 
0.04 dex, to be compared with the cr (, s = 0.07 (see Table [5]) 
we found for the same element in M 22. We find a rms in the 
calcium abundance almost egual to the one found in M 4 when 
we divide our stellar sample into the two s-process element rich 
and poor groups (<x i, s = 0.04 for both the two s groups). 

Apparently, the two stellar populations in M 4 and M 22 
have different origin, or the mechanisms responsible for this 
dichotomy must have been acting with different intensities in 
the two clusters. 

8. An independent check of the results 

In this paper we have presented a clear evidence of a spread 
in [Fe/H] and of the presence of a bimodal distribution in s- 
process element content among the M 22 stars. These results 
are based on high resolution UVES spectra of seventeen stars. 
Because of the high relevance of these results in the context of 
the ongoing lively debate on the multipopulation phenomenon 
in star clusters, we further searched in the ESO archive for addi- 
tional spectra of M 22 stars in order to strenghten the statistical 
significance of the results from UVES data. In fact, we found 
GIRAFFE spectra for 121 stars. We reduced all of them, but 
only fourteen were in the appropriate RGB location to have 
high enough S/N to pass all of our quality checks (see the 



following discussion) and useful for the abundance measure- 
ments. In any case we could double the original UVES sample 
of stars. 
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Fig. 18. ACS/WFC CMD of M 22 from Piotto (2009). 




Fig. 19. Red and blue lines are theoretical isochrones 
from Pietrinferni et al. (2004), with [Fe/H]=-1.68 and 
[Fe/H] =-1.82. The difference in metallicity produces on the 
CMD a split of the RGB and the SGB. 

Each star was observed with HR09, HR13, and HR15 
set-ups, which give a resolution of about R~20000-25000 in 
514-535, 612-640, and 660-696 nm respectivelly. Data were 
reduced using the last version of the pipeline developed by 
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Geneva observatory (Blecha et al. 2000), and were bias- 
subtracted, flat-field corrected, and extracted applying a wave- 
length calibrations obtained by ThAr lamps. Each spectrum 
was then normalized. Finally for each star we obtained radial 
velocity and applied a membership criterion as done for UVES 
data. Spectra of each member star were then reported to rest- 
frame velocity and combined together. We performed also a test 
in order to verify the influence of scattered light on the spectral 
line shape (and could alter the final abundances). To this aim we 
reduced some spectra covering the whole GIRAFFE CCD with 
and without scattered-light subtraction and compared spectral 
features. We found that scattered-light has no significant influ- 
ence on lines. 

For GIRAFFE data analysis, we wanted to follow the same 
procedure used for UVES, i.e. obtaining atmospheric parame- 
ters only from spectroscopy and not from photometry, which 
can be altered by differential reddening. This approach has the 
advantage to put the abundance determinations from the two 
data-sets in the same abundance scale, avoiding systematic ef- 
fects. However, this choice, coupled with the selection of stars 
located only in the RGB, was paid by the rejection of 103 spec- 
tra because of their low S/N(<70-80). Indeed an high S/N ratio 
is necessary to measure a sufficient number of isolated Fel/II 
lines in GIRAFFE spectra, which have lower resolution and 
cover a smaller wavelength range with respect to UVES ones. 
In addition, the brightest stars were not analized because of 
their very low temperature (T e ff < 4000). In this T s s regime 
also metal poor stars (as in the case of M 22) show very strong 
lines, which are blended in GIRAFFE data, not allowing a reli- 
able EW measurement. 

After a star passed through our selection criteria, atmo- 
spheric parameters (and Fe content) were obtained by EW 
method, as done for UVES and using the same line-list for the 
spectral lines in common between the two set-ups. 

Because of the two different data-sets, some comparisons 
are needed between the results obtained from the two spectro- 
graphs. Since we have only one star (#51) in common between 
the two data-sets, we could not compare directly (apart for this 
star) the results. In Tab. [7] we list the chemical abundances in 
common obtained for this star from the two different data. We 
note that the atmospheric parameters are in agreement within 
the errors calculated for UVES, the values for Fe and Y abun- 
dances are within the cr tot listed in Tab. |5J while for the other 
elements there are larger discrepancies probably due to errors 
associated to GIRAFFE results (that we have not considered 
here). Since the comparison of one star is not enough to verify 
the compatibility between the two set of abundances, we com- 
pared the atmospheric parameters and the mean metallicity. 

Figure [20] summarizes our tests. Filled circles represent 
GIRAFFE results, while open squares the UVES ones. In the 
upper left panel we plotted log(g) vs. T e ff , while the upper right 
panel shows v t vs. log(g). Both gravity as a function of tem- 
perature and microturbolence velocity as a function of gravity 
follow the same general trend for the two data-sets, with simi- 
lar dispersions. Further tests are shown in the two lower panels. 
The lower left one shows the Fe abundance vs. V magnitude 
(i.e. vs. the evolutionary state of the star along the RGB). No 
correlation appears, meaning that no systematic errors due to 



the different evolutionary phase are presents. The lower right 
panel reports T e ff vs. B - V color. The line is the empirical re- 
lation by Alonso et al. (1999), obtained assuming a reddening 
E(B - V)=0.34 (Harris 1996). Also in this case a good agree- 
ment was found, not only for the zero-point of the relations (i.e. 
the absolute average reddening of the cluster), but also for the 
shape of the relations themselves. 

An important test comes from the comparison of the mean 
iron content obtained from the two data-sets. From the fourteen 
GIRAFFE stars, we obtain: 

< [Fe/H] >= -1.74 + 0.03 

which perfectly agrees with the UVES value within 1 cr. 

A rough estimate of errors on atmospheric parameters can 
be done as in Marino et al. (2008), assuming that stars with the 
same V magnitude (corrected for differential reddening) have 
the same parameters. In this way we can obtain an upper (also 
the dispersion in metal content can in part contribute to the dis- 
persion of the stellar parameter values at a given luminosity 
along the RGB) estimate of the errors, which are AT e ff=±65 
K, Alog(g)=+0.20, and Av t =±0.11 km/s respectively. We can 
see that errors are a bit larger, but still comparable with UVES 
ones. 

All the other elements, with the exception of Ti, in 
GIRAFFE data were measured by spectral synthesis because 
of the severe blends with other lines. In addition to Fe and Ti 
we measured O (from the forbidden line at 630 nm), Na (from 
the doublet at 615 nm), Y (from the doublet at 520 nm), Ba 
(from the line at 614 nm), Nd (from the line at 532 nm), and Eu 
(from the line at 665 nm). 

Having verified the good agreement between the atmo- 
spheric parameters and the iron content obtained from the two 
data-sets, we can further proceed to verify whether GIRAFFE 
data confirms UVES results. For this reason, in Fig. |2T| we 
show some of the trends discussed in the previous sections. 
Filled circles are GIRAFFE measurements, while open circles 
are UVES ones. It is clear from this comparison that UVES re- 
sults are fully confirmed. In particular, we can confirm the Y-Ba 
bimodality (central panel), as well as the different Fe content 
for the two s-element rich and s-element poor groups of stars 
(see leftmost and rightmost middle panels). 

In addition, from GIRAFFE data we measured also Nd (a 
combined s and r element) and Eu (a pure r element) lines. 
Their abundances as a function of [Fe/H] are shown in the cen- 
tral and right lower panels. There is no trend for [Eu/Fe], while 
[Nd/Fe] clearly correlates with [Fe/H]. This is a further evi- 
dence that the iron enrichment of the s-process rich group is 
due to core-collapse SNe. 

9. Two stellar populations in M 22 

In the present paper we presented high resolution spectroscopic 
analysis for a sample of seventeen RGB stars in the GC M 22 
from UVES and FLAMES +UVES data. 

We confirm that M 22 is a metal poor GC, with a 
mean iron-content [Fe/H] =-1.75 + 0.02 (weighted mean be- 
tween UVES and GIRAFFE) and a mean a-enhancement 
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Fig. 20. Check of the atmospheric parameters used in this paper. Filled circles are GIRAFFE stars, while open squares are UVES 
ones. See Section 8 for a full explanation. 



[a/Fe]=+0.36 ± 0.04. Sodium and oxygen follow the well 
known anti-correlation, while no evidence for an Mg-Al an- 
ticorrelation was revealed. A clear correlation was found be- 
tween Na and Al. 

We find a strong dichotomy in the distribution of the s- 
process elements barium, yttrium and zirconium. Most impor- 
tantly, we find that the abundance of these elements is cor- 
related with the iron abundance. Stars enriched in s-process 
elements also show larger values of [Fe/H], by ~ 0.14 dex. 
The ^-process elements abundance correlates with calcium, and 
calcium with iron. No clear correlation is present between s- 
process elements and sodium, and between iron and sodium, 
but we noted that s-element and Fe enriched stars show higher 
values of sodium. These stars show also an overabundance of 
magnesium and silicon. The correlation between ^-process el- 
ements, and Ca abundances, with [Fe/H] is the strongest argu- 
ment in favor of the presence of two groups of stars with a dif- 
ferent Fe content in M 22. All these results have been confirmed 
by a sample of fourteen lower resolution GIRAFFE spectra, 
which allowed us to double the original UVES sample. 

According to the most recent theoretical models by 
Pietrinferni et al. (2004), a difference in metallicity of 0.14 
dex should cause a difference of ~0.10 mag in the F606W 
ACS/WFC band at the level of the SGB. Piotto (2009) have in- 



deed found that the SGB of M 22 is separated into two, distinct 
branches. However, the average separation in F606W band is 
of 0.17 magnitudes: it appears that the SGB split can not be 
attributed to a difference in [Fe/H] alone. 

The fraction of stars on the bright SGB (bSGB) corre- 
sponds to 62%±5% of the total SGB population, while the faint 
SGB (fSGB) includes the remaining 38%+5% of the SGB stars 
(Piotto 2009). In the stellar sample of the present paper, the 
fraction of Ba-strong, Y-strong, Zr-strong stars is ~41%. It is 
therefore tempting to connect the s-process element poor sam- 
ple to the bright SGB stars, while the faint SGB stars could 
be the ones with enhanced s-process elements. The correct re- 
production of the two SGBs needs an accurate determination 
of the NaCNO abundances, as shown by Cassisi et al. (2008) 
for the analogous case of NGC 1851. Also an He variations 
between the two populations can affect the SGB morphology. 
Indeed, we note that M 22 shares similarities with oj Centauri 
and NGC 1851: these clusters, where multiple stellar popu- 
lations have been photometrically identified along the SGB, 
exhibit a large range, not only in C, N, O, Na, Al, but also 
in i-process element abundance. M 22 is the only globular, 
apart lo Centauri, where some evidence of an intrinsic spread 
in iron were observed. Some hints (even if very uncertain due 
to the low number statistics of the analyzed sample) for a some 
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Fig. 21. Comparison of our final results we obtained for UVES (open circles) and GIRAFFE (filled circles). 



iron spread in NGC 1851 was suggested by Yong & Grundahl 

(2008) . NGC 1851, w Centauri, and M 22 show a large varia- 
tion in the Stromgren index, traditionally used as a metallicity 
indicator. All of these three GCs have a splitted SGB. 

From our observations and from the results of Piotto 

(2009) , it is tempting to speculate that the SGB split could 
be related to the presence of two groups of stars with differ- 
ent s-process element content and a difference, albeit small, in 
iron. According with this scenario, s-process element poor stars 
are those populating brighter SGB stars and constitute the first 
M 22 population. The second stellar generation should have 
been formed after that the AGB winds of this first stellar gen- 
eration have polluted the protocluster interstellar medium with 
i-process elements. This second generation may have formed 
from material which was also enriched by core-collapse super- 
novae ejecta, as indicated by the higher iron, magnesium, and 
silicon content, and the lack of correlation of the iron content 
with a pure r-process element (Eu). A detailed analysis of the C, 
N, O abundances of the SGB stars in M 22 is strongly needed 
in order to properly settle this problem. 
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